Specific heat of Bosons in a lattice 
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We present a theoretical study ol specific heat ol bosons (C„) in a simple cubic lattice. We 
have studied the non-interacting bosons and the Tonks gas. For both cases, the Cv above the bose 
condensation temperature shows considerable temperature dependence compared to that of free 
bosons. For Tonks gas, we find that the low-temperature specific heat increases as the system gets 
closer to the Mott transition. 
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I. INTRODUCTION 

Investigations of collective properties of bosons in opti- 
cal lattices produced by counter propagating laser beams 
is now a rapidly expanding field of research. It was 
shown^ early on that bosons in optical lattices can be ad- 
equately modeled by employing a Bose-Hubbard model. 
Properties of a system modeled by this model depends 
on lattice symmetry, inter-site boson hopping energy 
(t), on-site boson-boson interaction strength (U), and 
the number of bosons per site (n). Experimentalists 
have achieved great control over several of these param- 
eters by changing the characteristic of the laser beams 
and their relative orientations. Recently, Greiner and 
collaborators^ demonstrated bose-condensed to Mott in- 
sulator transition in a system of interacting bosons in a 
simple cubic optical lattice. This transition was predicted 
in theoretical studies^'^"^ on Bose-Hubbard model. 

For a system of repulsive interacting bosons living in 
a lattice, bosons will try to avoid multiple occupancy of 
the sites. When the number of bosons per site is less 
than or equal to unity and in the strong interaction limit 
{U » t) double (or multiple) occupancy of any site 
is forbidden. In this case, one has a quantum version 
of Tonks gas* (here after referred to as Tonks gas) in 
a lattice. One dimensional Tonks gas was investigated 
by Girardeau^ long time back. Quite recently Tonks- 
Girardeau gas^'^'^ regime has been achieved in experi- 
mental studies"'^^ of interacting bosons in one-dimensional 
optical lattices. It is reasonable to expect that Tonks gas 
in two and three dimensional lattice will be achieved in 
the near future. Thermodynamic properties of bosons in 
a lattice is of considerable interest in this context. In 
this short communication we present a theoretical study 
of specific heat of non-interacting and strongly interact- 
ing bosons in a simple cubic (sc) lattice. As should be 
clear by now, the experiments mentioned above and ex- 
periments that would be conducted in future are our pri- 
mary motivation for the work presented in this paper. 
Heat capacity measurements on an optically trapped, 
strongly interactiong Fermi gas of atoms has been re- 
ported recently^^ and one can expect that such measure- 
ments for bosons in an optical lattice will be done in the 



near future. The effects of the presence of the harmonic 
confining potential is not included in our calculations 
given below. Our results for the non-interacting bosons 
would be reasonable for a trap of shallow confining po- 
tential. It may be mentioned, however, that the strongly 
interacting regime (Tonks limit) could be reached in ex- 
periments with extremely deep optical lattice potential. 
This, in turn, creates a deep harmonic trap, generated by 
the envelope of the laser beams. For realistic calculations 
of experimentally measured quantities in the Tonks limit 
at present, one should take account of the harmonic trap 
potential. Our study in the Tonks limit may be consid- 
ered as a starting point of a more complete calculation 
where the effects of inhomogeneties produced by the con- 
fining potenetial is included. 



II. SPECIFIC HEAT OF LATTICE BOSONS 

The specific heat can be obtained from the energy 
given by 



E = 



w 



(1) 
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where p(e) is the normalized single boson Density 
of States (DOS) of a band of width 2W, ngie) = 
l/[exp{P{e — /i) — 1] with f3 = l/fc^T where kg is the 
Boltzmann's constant and T is the temperature, and fi 
is the chemical potential. 

First we consider the case oi T > Tb , where Tb is the 
Bose condensation temperature. In this case, the chem- 
ical potential is below the band bottom and is temper- 
ature dependent. The chemical potential and the bose 
condensation temperature are obtained by solving the 
number equation 
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p{e)nB{e)de, 
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where n is the number of bosons per site. For T > Tb, 
the energy is given by Eq. (1) and we obtain the constant 
volume specific heat (per site) to be 



(3) 



where, 

lo = F{e- n)p{e)de, (4) 
J-w 

/w 
eF{e - fi)p{e)de, (5) 
-w 

/w 
{e-^i)F{e-^i)p{e)de, (6) 
-w 

h= e{e-^l)F{e-^l)pie)de. (7) 
J-w 

In the above equations, 

F{e-f,) 



<('-/') 



[e/3(e-/.) _ 1] 
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Now we consider the case of T < Tg. In this case 
the chemical potential is temperature independent and 
is pinned to the bottom of the band {—W). The band 
bottom has a macroscopic number of particles (the Bose 
condensate) and the rest of the particles are in the ex- 
cited states. Denoting next as the number of bosons (per 
site) in excited states, the energy per site is 

E= ep{e)nB{e)de-W{n-next)- (9) 
J-w 

The specific heat in this case becomes, 
1 

^- = k^ + wfp^^)P(^ + (10) 

Non-interacting bosons: The Hamiltonian of the non- 
interacting bosons in a tight-binding energy band is: 



F = 5^[e(k)-M]cU. 



(11) 



where, esc{kx,ky,kz) = -2t X^^^^ cos(fc^) is the energy 
band structure of a single boson in a s c lattice with lat- 
tice constant sot to unity, Ck ~ y^l/Ns exp{i'k.'Ri) , 
Ci is the annihilation operator of a boson at site i, and 
Ns is the number of lattice sites. 

Strongly interacting bosons (the Tonks gas limit): To 
study the Tonks gas limit, we start from the Bose- 
Hubbard model given by 



H = -t^ {clcj + h.c.) + — ^rtiirii - 1) - p^rii , 

<ij> i i 

(12) 



where rii = clci. In the above Hamiltonian c? (cj) is 
the boson creation (annihilation) operator, U the boson- 
boson repulsive interaction energy. The Tonks gas regime 



is achieved in the strong correlation limit {U/t » 1) 
when the number of particles per site n = N/Ns is re- 
stricted to be less than or equal to unity. In this limit, 
the double (or multiple) occupancies of the sites are for- 
bidden. We will follow a Rcnormalizcd Hamiltonian Ap- 
proach (RHA) similar to that used for strongly interact- 
ing fermions^^'^^ based on Gutzwiller approximation^^. 
In this method, the strict constraint of no double occu- 
pancy shows up in the rcnormalization factor for hopping. 
The hopping (or band) rcnormalization factor is obtained 
by taking the ratio of hopping probability between two 
sites in the correlated space to that in the non-correlated 
space. The hopping probability in the projected space 
(of no multiple occupancy) is 



Pc 



' n{l — n). 



(13) 



The above equation means that the site from which hop- 
ping takes place should be occupied and the target site 
should be empty in the projected space. To find the hop- 
ping probability in the non-correlated state, we have to 
first find the probability that a site is occupied. The oc- 
cupancy of the site to which a boson hops is irrelevant 
since we are dealing with bosons. To find the site oc- 
cupancy probability of a site from which a boson hops, 
we employ the following route. First we find the num- 
ber of ways Ng number of bosons can be distributed in 
Ns number of lattice sites and then the number of ways 
of distributing the same number of bosons in {Ns — 1) 
number of lattice sites. The difference would give the 
number of configurations where a particular site among 
Ns number of lattice sites is occupied by bosons. The 
number of configurations possible for Nb bosons on Ns 
sites is given by 



w{Nb : Ns) 



{Nb + Ns- 1)! 
Nb\{Ns- 1)1 ■ 



(14) 



The number of configurations possible for Nb bosons on 
{Ns — 1) sites is given by 

tAT AT ^^ (iVg + Ns- 2)! 
w{Nb :Ns-1)= ^^,(^^_,), ■ (15) 

The probability that a particular site is occupied is then 



_ ^ w{Nb :Ns-1) 
w{Nb:Ns) • 



In the thermodynamic limit. 



Po 



1 + n 



(16) 
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As mentioned for non-interacting bosons, the hopping 
takes place irrespective of the target site being empty or 
occupied by any number of bosons. Hence the hopping 
probability for the non-interacting case is 



Pn 



1 + n 



(18) 



The band-width renormalization factor {Pcorr / Pnocorr) 
then is, 



= l-n'. 



(19) 



Using the above cf)B the renormahzed Hamiltonian vaUd 
in the strong couphng hmit and for n < 1 is, 



Hsc = ^[(^s(n)e(k) - M]4ck- 



(20) 



The above is the Hamiltonian of Tonks gas on a lat- 
tice. One can immediately see that the band width is 
strongly n dependent. For n = 1, the effective mass 
diverges and the system is a Bose-Mott-Hubbard insu- 
lator. In passing , we note that for the fermion Hub- 
bard model, the corresponding hopping renormalization 
factor^^ jg 2(1 _ n)/{2 - n). We do realize that the RHA 
approach based on Gutzwiller approximation used to ob- 
tain Eq. 20 has its limitations. Within this approach, 
the effect of strong interaction shows up only in the en- 
hanced effective mass of the bosons. This approxima- 
tion would be reasonable for large connectivity (i.e., for 
large number of nearest neighbors). One can get some 
feel for the accuracy of the approximation by compar- 
ing with the same for fermions. In the case fermions, it 
has been shown that the Gutzwiller approximation gets 
progressively better as the dimensionality increases and 
becomes exact in infinite dimensions and that in three 
dimensions it is reasonable^^. For strongly interacting 
bosons also, the Gutzwiller wave function is believed to 
be exact for large dimensionality^'' and can be expected 
to be reasonable in three dimensions. 

As mentioned earlier the chemical potential and the 
Bose condensation temperature (Tb) is determined by 
solving the number equation. Results of this calculations, 
for bosons in a sc lattice, are shown in Fig. 1. For the 
non- interacting bosons, the Tb shows an initial fast rise 
and then almost a linear increase with n. In the case of 
strongly interacting bosons, correlation-effects lead to an 
enhancement of the boson effective mass with increasing 
band-filling. This increase in the effective mass becomes 
dominant beyond n = 0.5. As a consequence the Tb 
starts decreasing and finally going to zero at integer fill- 
ing for which the transition to Bose-Mott-Hubbard insu- 
lating state takes place. It may be noted that the Tb for 
strongly interacting bosons is almost the same as that for 
non-interacting bosons for low filling (n < 0.2). This is a 
signature of the fact that correlation effects are not signif- 
icant for low density. The specific heat is calculated using 
the Eqs. (3) and (10). For the Tonks gas calculations, 
the original band-width is replaced by the renormalized 
band-width. The results for non-interacting bosons in a 
sc lattice and a comparison with that for strongly inter- 
acting bosons are shown in Fig. 2. In general, the spe- 
cific heat increases with temperature for T <Tb, shows 



a A anomaly at Tb and then decreases with temperature. 
For non-interacting bosons, the specific heat (per site) 
follows the same curve at low temperatures irrespective 
of the band filling. The peak value of the specific heat, 
however, depends on the filling. The peak value is larger 
for higher filling. It is also found that the peak value 
does not depend on the bandwidth. For strongly inter- 
acting bosons Fig. 2 shows that for a fixed boson density 
(i) the peak in the specific heat occurs at a temperature 
lower than that for the non-interacting bosons, (ii) the 
peak value, however, remains the same for both inter- 
acting and non- interacting cases. The latter result is a 
consequence of the fact that the peak value of the specific 
heat depends only on the boson density and not on the 
boson effective mass (or bandwidth). The first behavior 
is a combined effect of the reduced value of Tb for in- 
teracting bosons compared to the non-interacting bosons 
at the same density and the appearance of the specific 
heat peak at Tb- As mentioned previously the effect 
of the on-site interation is small for low boson density 
and becomes increasingly dominant as the boson density 
approaches to integer filling, the difference between the 
specific heat curves as well as between the peak temper- 
atures for the interating and non-interacting bosons is 
small for n = 0.4 and is large for n = 0.8 where the 
correlation effect is dominant. In experiments it is now 
possible to tune the hopping of bosons for a given density 
of bosons to control the effective interaction strength and 
to obtain results for weakly (or non-interacting) bosons 
and strongly interacting bosons. The ratio of the peak 
temperatures of the specific heat for the non-interacting 
and interacting cases for a fixed boson density may serve 
as a measure of effective interaction. A comparison of 
specific heats for non-interacting bosons in a sc lattice 
and that for bosons in a 3-d box is shown in Fig. 3. 
The specific heat for free bosons in a 3-d box is obtained 
by using the DOS p{e) oc -/e. For the sc case the spe- 
cific heat follows a S-shapcd curve for T < Tb and it 
shows a considerable temperature dependence above Tb 
compared to that of free bosons. The A anomaly is very 
prominent for the sc case in comparison with the free 
bosons. The results for strongly interacting bosons are 
shown in Fig. 4. For the strongly interacting bosons, 
it is seen that with increasing n the specific heat versus 
temperature curves become narrower. At a low temper- 
ature the specific heat as well as the slope of the specific 
heat with temperature increases with increasing n. This 
is because of increasing effective mass of the bosons with 
n. 



III. CONCLUSIONS 



In this paper we have presented specific heat for bosons 
in a tight-binding band corresponding to a sc lattice. We 
investigated non-interacting bosons and Tonks gas. In 
both cases, the specific heat is found to have consider- 



ably more temperature dependence above Tb compared 
to that of free bosons in a box. For the strongly inter- 
acting bosons, we find that the low-temperature specific 
heat gets enhanced as the systems moves towards the 
Mott transition. 
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FIG. 1. Bose condensation temperature for non-interacting 
bosons (dots) and strongly interacting bosons (solid line) as 
a function of number density (n) for bosons in a sc lattice. 
Here W is the half-band-width. 
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FIG. 2. Temperature variation of specific heat per site for 
non-interacting bosons (thick lines) in a sc lattice for two val- 
ues of n (0.4 and 0.8). Thin lines are for strongly interacting 
bosons, (a) n = 0.8 and (b) n = 0.4. 



